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ABSTRACT 

Gaseous and stellar metallicities in galaxies are nowadays routinely used to constrain 
the evolutionary processes in galaxies. This requires the knowledge of the average yield 
per stellar generation, yz, i.e. the quantity of metals that a stellar population releases 
into the interstellar medium (ISM), which is generally assumed to be a fixed fiducial 
value. Deviations of the observed metallicity from the expected value of yz are used to 
quantify the effect of outflows or inflows of gas, or even as evidence for biased metal- 
licity calibrations or inaccurate metallicity diagnostics. Here we show that yz depends 
significantly on the Initial Mass Function (IMF), varying by up to a factor larger than 
three, for the range of IMFs typically adopted in various studies. Varying the upper 
mass cutoff of the IMF implies a further variation of yz by an additional factor that 
can be larger than two. These effects, along with the variation of the gas mass fraction 
restored into the ISM by supernovae [R, which also depends on the IMF), may yield 
to deceiving results, if not properly taken into account. In particular, metallicities that 
are often considered unusually high can actually be explained in terms of yield asso¬ 
ciated with commonly adopted IMFs such as the Kroupa (2001) or Chabrier (2003). 
We provide our results for two different sets of stellar yields (both affected by specihc 
limitations) finding that the uncertainty introduced by this assumption can be as large 
as ^ 0.2 dex. Finally, we show that yz is not substantially affected by the initial stellar 
metallicity as long as Z > 10“^ Zq. 

Key words: galaxies: evolution - galaxies: ISM - ISM: abundances - stars: abun¬ 
dances 


1 INTRODUCTION 

The analysis of the chemical enrichment of galaxies is a pow¬ 
erful tool to constrain galaxy evolutionary processes. The 
content of metals in galaxies, both in the interstellar medium 
(ISM) and in stars, depends critically on the past star for¬ 
mation history and on the net effect of outflows and inflows, 
which are some of the key mechanisms in shaping galaxy 
evolution. In order to extract valuable information from the 
observed metallicities, it is crucial to compare them with 
the amount of metals expected to be produced by the inte¬ 
grated star formation. To achieve this is necessary to have 
accurate information on the amount of each chemical ele¬ 
ment injected into the ISM by each type of star, i.e. the 
so-called stellar yields. Generally, most observations provide 
information only on the global content of metals, or on a 
single chemical element which is taken as representative of 
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the global metallicity (by assuming that, for instance, the 
abundance of the various elements scales proportionally to 
the solar relative abundances). Moreover, many studies do 
not deal with the relative delayed enrichment of different 
chemical species. Therefore, the quantity that is often used 
is the so-called average yield per stellar generation, or net 
yield (generally indicated as yz, or simply y), which is the 
total mass of metals that a stellar population releases into 
the ISM, normalized to the mass locked up into low-mass 
(long-lived) stars and stellar remnants. 

Historically, this approach was first used in the early 
work of Searle & Sargent (1972), who derived the relation 
between gas phase metallicity Z and gas fraction = 
Mgas /(Mgas + M*) for a closed box model (see Tinsley 1980, 
or Matteucci 2001 for a detailed analysis); 



This simple model is based upon the assumptions that the 
galaxy is one-zone and closed; the initial gas mass is of 
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Figure 1. In this figure, we report the ejected mass of oxygen 
as a function of the progenitor mass, for different initial stel¬ 
lar metallicities. This set of stellar yields is from Romano et al. 
(2010). The solid line in black corresponds to the stellar yields at 
Z = 1.0 X 10“^°; the dashed line in red to Z = 1.0 X 10“®; the 
black dotted line to Z = 1.0 X 10“®; the dashed-dotted line in 
blue to Z = 3 X 10“®, and the green squares to the stellar yields 
at Z = 2.0 X 10“2. 


primordial chemical composition; the initial mass function 
(IMF) is invariant, and the mixing of the gas in the galaxy 
is always instantaneous and complete, and that metals are 
instantaneously recycled for the formation of the new gener¬ 
ation of stars. The latter is dubbed “instantaneous recycling 
approximation” (IRA, Tinsley 1980). Within this simplified 
(but widely used) approach the further approximation is 
that all stars with m > 1 Mq die instantaneously, while 
all stars with m < 1 Mq have infinite lifetime; this way, the 
effect of stellar lifetimes in the equations can be neglected 
and the effect subsumed into a net return fraction {R). 

Although the IRA assumption is strong, it still rep¬ 
resents a good approximation for those chemical elements 
produced and restored into the ISM by stars with short 
lifetimes. The best example of such a chemical element is 
oxygen, which is also representative of the global metallic- 
ity Z, since it is the most abundant heavy element by mass. 
On the other hand, the ISM evolution of chemical elements 
produced by long-lifetime sources cannot be followed by ana¬ 
lytical models working with the IRA assumption. Examples 
of such chemical elements are nitrogen and carbon, which 
are mainly synthesized by low- and intermediate-mass stars 
(LIMS), and iron, mainly produced by Type la SNe. To take 
into account the stellar lifetimes with a high level of detail, 
numerical models of chemical evolution should be used (see 
Matteucci 2012). 

The yield per stellar generation is a key in the con¬ 
text of analytical models of chemical evolution, even in 
the most complex ones, which include the effect of out¬ 
flows and inflows, as well as variations of the star formation 
efficiency, i.e. normalization and slope of the relation be¬ 
tween gas mass and star formation rate, (see, for example, 
Bouche et al. 2010; Spitoni et al. 2010; Dekel et al. 2013; 
Peng & Maiolino 2014; Zahid et al. 2014; Ascasibar et al. 
2014; Spitoni 2015), as well as in numerical simulations (e.g. 



Figure 2. In this figure, we show how ejected mass of oxygen 
vary as a function of the progenitor mass, for different initial stel¬ 
lar metallicities, when assuming the Nomoto et al. (2013) stellar 
yields. The solid line in black corresponds to the stellar yields at 
Z = 0.0; the dotted line in black to Z = 1.0 X 10“®; the dashed- 
dotted line in blue to Z = 4.0 X 10“®; the solid line in magenta 
to Z = 8.0 X 10“®; the solid line in green to Z = 2.0 X 10“^, and 
the solid line in cyan to the stellar yields at Z = 5.0 X 10“^. 


Cole et al. 2000; De Lucia et al. 2004). The comparison of 
these models with the extensive observations that are pro¬ 
viding metallicity measurements for large samples of galax¬ 
ies locally and at high redshift (e.g. Savaglio et al. 2005; 
Maiolino et al. 2008; Troncoso et al. 2014; Steidel et al. 
2014) enable us to provide important constraints on these 
various processes, modulo an accurate knowledge of the yield 
per stellar generation. 

In most of the studies, the yield per stellar generation 
is taken as a fixed value (typically about 0.012-0.045), with 
this value changing from work to work. For example, a net 
yield of oxygen yo = 0.015 is assumed in Peeples et al. 
(2014), whereas yo = 5.7 x 10~® in Zahid et al. (2014) 
and yo = 3.13 x 10“® in Ho et al. (2015); particularly high 
are the values yo = 0.040-0.045 assumed in De Lucia et al. 
(2004), yo = 0.03 in Croton et al. (2006), and yo = 0.04 in 
Bower et al. (2008). However, since the net yield is a combi¬ 
nation of yields from different stellar masses, it is clear that 
it must depend on the IMF. This has sometimes been ac¬ 
knowledged (e.g. Henry et al. 2000; Kobayashi et al. 2011; 
Ho et al. 2015), but never really taken into consideration 
when using the net yield in the various models. In particu¬ 
lar, several works derive the stellar mass and SFR by assum¬ 
ing a given IMF and then adopt a fiducial net yield that is 
derived from a completely different IMF. Moreover, there is 
some evidence that the IMF may vary in different classes of 
galaxies. This implies that different yields per stellar genera¬ 
tion should be used. Finally, since the stellar nucleosynthetic 
yield have a metallicity dependence, it is important to check 
the effect of metallicity on the IMF-integrated net yield. 

To tackle the issues presented above, in this article we 
calculate yields per stellar generation for the most commonly 
adopted IMFs and investigate their metallicity dependence 
(although the latter effect is shown to be minor), by com¬ 
paring the results for two modern compilations of nucle- 


MNRAS 000, 1-9 (2015) 















Modern yields per stellar generation 3 


osynthetic yields (Romano et al. 2010; Nomoto et al. 2013), 
which have been throughly tested in the past against the 
best available data for galaxies, although we remark on the 
fact that each of them is still affected by specific limitations 
(as discussed in the following Sections). We mostly focus 
on the yield of oxygen, which is the element which is most 
commonly used as a tracer of the global metallicity, and for 
which the IRA approximation is appropriate. However, we 
will also provide the yield per stellar generation for the total 
mass of metals, although this should be used with caution, 
given the enrichment delay of various elements (e.g. iron, 
nitrogen, carbon, etc...), for which the IRA approximation 
is arguable. 

In Section 2 we define the quantities we have computed 
in this work and specify the set of stellar yields which we 
have assumed and the IMFs which we have explored. In Sec¬ 
tion 3 we report and discuss our results. Finally, in Section 
4 we summarize the main conclusions. 


2 DEFINITIONS AND ASSUMPTIONS 

2.1 Yield per stellar generation and retnrn mass 
fraction 


We define as yield per stellar generation i/i(Z), or net yield, 
the ratio of the global gas mass in the form of a given chem¬ 
ical element i newly produced and restored into the ISM 
by a simple stellar population with initial metallicity Z to 
the amount of mass locked up in low mass stars and stellar 
remnants (Tinsley 1980; Maeder 1992; Matteucci 2001): 

1 V mPi{m,Z)<j)(m)dm 

= 1 -> (2) 

where: 

(i) pi(m,Z) = is so-called stellar yield, 

which is defined such that m ■ Pi(m, Z) represents the mass 
in the form of the i-th chemical element newly formed and 
ejected into the ISM by stars with initial mass m and metal¬ 
licity Z; 

(ii) ((>(m) is the IMF, namely the mass-spectrum over 
which the stars of each single stellar generation are dis¬ 
tributed at their birth; 

(iii) miong-iiv = 1.0 Mq is the maximum mass of the so- 
called long-lived stars, which do not pollute the ISM; 

(iv) rriup is the upper mass cutoff of the IMF; in our stan¬ 
dard case, we assume mup = 100 Mq , however in the second 
part of the article we will also investigate the effect of vary¬ 
ing m-up. 


Finally, R represents the so-called return mass fraction, 
which is defined as the total mass fraction (including both 
processed and unprocessed material) returned into the ISM 
by a stellar generation: 


R(Z) 


^ng-liv ~ 

/on“MQ 


( 3 ) 


with Z) being the mass of the stellar remnant left by 

a stars with initial mass m and metallicity Z. 

If one changes the quantity miong-iiv in accordance to 
the age of the galaxy, then one would obtain equations very 
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Figure 3. In this figure, we show the trend of the different 
IMFs studied in this work, as normalized with respect to the 
Salpeter (1955) IMF. The dotted line in black corresponds to the 
Kroupa et al. (1993) IMF; the dotted-dashed line in red to the 
Kroupa (2001) IMF; the dashed line in blue to the Chabrier (2003) 
IMF, and the solid line in purple corresponds to the Salpeter 
(1955) IMF. 
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Figure 4. In this figure, we compare the Cappellari et al. 
(2012, blue dashed line) top-heavy IMF with the IMFs of 
Arimoto & Yoshii (1987, red dotted line) and Salpeter (1955, solid 
line in magenta). As in Fig. 3, all the IMFs are normalized with 
respect to the Salpeter (1955) IMF. 


similar to the ones numerically solved by current models 
of chemical evolution (see, for example, Matteucci 2012). 
In principle, the assumption of miong-Uv = 1.0 Mq provides 
correct results only for stellar populations which are older 
than « 7.1 Gyr, which corresponds to the lifetime of an 1 Mq 
star, according to Padovani & Matteucci (1993), although 
the lifetimes of low-mass stars can also be influenced by 
metallicity, particularly at very low Z (Gibson 1997). 
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Table 1. In this Table, we report the numerical values of (Po)i which is defined as the IMF-averaged stellar yield of oxygen in the mass 
range M = 10-40Mq, for different metallicities Z. ’MM02’ and ’HMM05’ stand for Meynet & Maeder (2002) and Hirschi et al. (2005), 
respectively, which correspond to the results of the Geneva stellar models also included in Romano et al. (2010). ’NKT13’ stands for 
Nomoto et al. (2013). 


Stellar yields 

Vrot [km s 1] 

z 

(Ho) 

(Ho) 

(Ho) 

(Ho) 




Salpeter (1955) 

Chabrier (2003) 

Kroupa et al. (1993) 

Kroupa (2001) 

MM02 

0 

4.0 X 10-3 

0.004 

0.007 

0.002 

0.007 


300 


0.006 

0.010 

0.003 

0.010 

HMM05 

0 

2.0 X 10-2 

0.006 

0.010 

0.003 

0.010 


300 


0.009 

0.015 

0.005 

0.014 

NKT13 

no 

4.0 X 10-3 

0.007 

0.011 

0.004 

0.010 


no 

2.0 X 10-2 

0.006 

0.009 

0.003 

0.009 


2.2 Stellar yields and initial mass function 

We use a numerical code of chemical evolution to explore 
the effect of the metallicity and IMF on the yields of oxygen 
per stellar generation and on the return mass fraction. We 
provide our results for the following sets of stellar yields: 

(i) the set provided by Romano et al. (2010), which as¬ 
sume the stellar yields of Karakas (2010) for LIM stars, 
and the He, C, N and O stellar yields of the Geneva 
stellar models for massive stars (Meynet & Maeder 2002; 
Hirschi et al. 2005; Hirschi 2007; Ekstrom et al. 2008); for 
heavier elements, which are not relevant for this study, 
Romano et al. (2010) assume the stellar yields of massive 
stars of Kobayashi et al. (2006); 

(ii) the stellar yields of Nomoto et al. (2013), which in¬ 
clude the stellar yields of LIM stars of Karakas (2010), and 
the stellar yields of Nomoto et al. (2006), Kobayashi et al. 
(2006), Kobayashi et al. (2011), Tominaga et al. (unpub¬ 
lished) for core-collapse supernovae (SNe). 

The mass of the stellar remnants have been collected by 
Romano et al. (2010) from the work of Kobayashi et al. 
(2006). Nevertheless, according to the Talbot & Arnett 
(1973) formalism, we compute the return mass fraction with 
the Romano et al. (2010) stellar yields, by summing the 
ejecta of all the chemical elements (both the processed and 
the unprocessed ones) for each stellar mass, and this quan¬ 
tity turns out to be dominated by the H and He contribu¬ 
tions. 

The Romano et al. (2010) compilation of stellar yields 
for He, C, N and O include the results of models which take 
into account the combined effect of mass loss and rotation 
(see also Maeder 2009 for a detailed discussion), whereas 
the Nomoto et al. (2013) stellar yields have been computed 
by models which do include standard mass loss but not the 
effect of rotation. With standard mass loss, only C and N 
have been lost before supernova explosions. However, the 
mass loss driven by rotation turns out to be particularly im¬ 
portant at almost solar metallicity and above in depressing 
the oxygen stellar yields of the most massive stars (M > 30- 
40 Mq, see also Fig. 1). In fact, mass loss increases with 
stellar metallicity and stars of high metal content loose H, 
He, but also C, through radiatively line driven winds. There¬ 
fore, the C production is increased by mass loss whereas the 
oxygen production is decreased, since part of C which would 


have been transformed into O, is lost from the star (Maeder 
1992). Finally, the effect of rotation is to produce mixing and 
enhances mass loss, with the efficiency of the mixing process 
being larger at lower metallicities (see also Chiappini et al. 
2008, and references therein). 

We remark on the fact that Romano et al. (2010) com¬ 
bine results of stellar models assuming only hydrostatic 
burning and rotation (Geneva group, for He, C, N, and O) 
with the results of models including explosive burning with¬ 
out rotation (Nomoto group, for heavier elements), giving 
rise to an inhomogeneous set of stellar yields, which can be 
physically incorrect. In the context of this study, the treat¬ 
ment of Romano et al. (2010) has a marginal effect, since the 
metallicity is dominated by the oxygen and carbon contribu¬ 
tions. On the other hand, Nomoto et al. (2013) provide one 
of the most homogeneous set of stellar yields available at the 
present time, although it is still affected by the limitation of 
not including the effect of stellar rotation. 

In Figures 1 and 2, we show how the oxygen stellar 
yields of Romano et al. (2010) and Nomoto et al. (2013), re¬ 
spectively, vary as functions of the initial stellar mass and 
for different metallicities. The Geneva stellar yields are avail¬ 
able only up 60 Mq and we therefore assume in our stan¬ 
dard case that the yields from 60 to 100 Mq are constant. 
On the other hand, the stellar yields of massive stars of 
Kobayashi et al. (2006, included in Romano et al. 2010, for 
the elements heavier than oxygen) and Nomoto et al. (2013) 
are available only up to 40 Mq and thus we keep them con¬ 
stant for stars with larger initial mass. We remark on the 
fact that very massive stars are expected to leave a black 
hole as a remnant; therefore, a significant fraction of stellar 
nucleosythetic products in the ejecta of very massive stars 
may eventually fall back onto the black hole. This process 
might cause a reduction of the stellar yields of very massive 
stars. 

At solar metallicity, Romano et al. (2010) include stel¬ 
lar yields which have been computed by applying a stellar 
rotational velocity Urot = 300 km s“^. From an observa¬ 
tional point of view, Ramirez-Agudelo et al. (2013) found 
that almost 80 per cent of nearby stars rotate slower than 
Urot = 300 km s~^, which thus can be considered as an 
approximate upper limit. To quantify the effect of stellar ro¬ 
tation in the stellar yields of oxygen from massive stars, in 
Table 1 we compare the predictions of models with and with- 
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out stellar rotation, with the quantity (Po) being defined as 
the IMF-averaged yield of oxygen in the mass range M = 10- 
40 Mq . The effect of stellar rotation in the Geneva stellar 
models is to increase the average oxygen stellar yield by a 
factor of ~ 1.5 for stars with initial mass below 40 Mq (see 
also Hirschi et al. 2005). Furthermore, at Z = 4.0x 10“®, the 
IMF-averaged oxygen stellar yield of Nomoto et al. (2013) - 
which neglect the effect of stellar rotation - is larger than 
the value of the Geneva stellar models without stellar ro¬ 
tation, but rather similar to the corresponding value with 
rotation; conversely, at solar metallicity, the Geneva stellar 
models without rotation agree with Nomoto et al. (2013). 
IMFs containing a larger number of massive stars, such as 
the Ghabrier (2003) and Kroupa (2001) ones, amplify the 
oxygen enrichment of the ISM and give rise to larger values 
of (Po), whatever be the set of stellar yields assumed. 

In this article, we study the effect of different IMFs: 
the Salpeter (1955), the Kroupa et al. (1993), the Kroupa 
(2001), and the Ghabrier (2003) IMFs, which are shown in 
Fig. 3 as normalized with respect to the Salpeter (1955) 
IMF. We have chosen these IMFs since they have been the 
most widely used by various authors. Moreover, these IMFs 
give quite different weights to different stellar mass ranges, 
hence they will more clearly display differences in the final 
predicted net yields and return mass fractions. As one can 
notice from Fig. 3, the Kroupa et al. (1993) IMF predicts 
the largest fraction of intermediate mass stars, while hav¬ 
ing the lowest number of massive stars. On the other hand, 
the Ghabrier (2003) and the Kroupa (2001) IMFs predict a 
higher number of both intermediate-mass stars and massive 
stars than the Salpeter (1955) IMF. 

Finally, we also explore the effect of the top-heavy IMFs 
of Gappellari et al. (2012) and Arimoto & Yoshii (1987), 
which are shown in Fig. 4, as normalized with respect 
to the Salpeter (1955) IMF. These IMFs are defined as 
a single-slope power law: <l>{m) oc with the 

Gappellari et al. (2012) IMF having a slope x = 0.5 and 
the Arimoto & Yoshii (1987) IMF assuming x = 0.95. 


3 RESULTS 

In this Section we present the net nucleosynthetic yields and 
return fractions obtained with the IMFs and stellar yields 
discussed in the previous section. As mentioned above, we 
will focus on the yield of oxygen, since it is the element 
most commonly measured and taken as representative of 
the bulk of the metallicity, and also because it is an element 
for which the IRA approximation is appropriate. However, 
we will provide a value also for the yield of the total mass 
of metals, although with some cautionary warnings. 

In Table 2, we show how the net yield of oxygen per 
stellar generation varies as a function of the IMF and metal¬ 
licity. In our “fiducial” case, reported in Table 2, we assume 
m-up = 100 Mq. Concerning the dependence on metallicity, 
the most interesting result is that the yield yo is roughly 
constant down to very low metallicities. This result implies 
that the assumption of a time-independent net oxygen yield, 
as generally treated in analytical models, is a reasonable 
one. Interestingly, we find an enhancement of yo for metal- 
free stellar populations (case with Z = 1.0 x 10~^° from 
Ekstrom et al. 2008). In fact, it is well established that, at 



mup/Me 

Figure 5. In this figure, we show how {yo) is predicted to vary 
as a function of rriup, which is defined as the upper mass cutoff of 
the IMF. We have computed {yo)-> by averaging the net yield of 
oxygen over the metallicity range 1.0 x 10“^ ^ Z < 2.0 x 10“^, 
within which yo turns out to be nearly constant. Thick lines 
represent our results for the Romano et al. (2010) stellar yields, 
while thin lines represent the Nomoto et al. (2013) stellar yields. 
The various curves with different colours correspond to the same 
IMFs as in Fig. 3. 



rriup/Me 

Figure 6. In this figure, we report how the average net yield of 
oxygen, {yo)', varies as a function of the upper cutoff of the IMF, 
when assuming the top-heavy IMF of Gappellari et al. (2012, blue 
dashed line) and the IMFs of Arimoto & Yoshii (1987, red dotted 
line) and Salpeter (1955, solid line in magenta). As in Fig. 5, {yo) 
is computed by averaging yo over the metallicity range 1.0 x 
10“^ ^ Z < 2.0 X 10”^, and the various curves with different 
thickness correspond to the same stellar yields as in Fig. 5. 


very low Z, the mixing induced by rotation is particularly 
efficient (Chiappini et al. 2008); in this way, the nucleosyn¬ 
thetic products of the 3a reaction in the inner He-burning 
zone can diffuse to the outer stellar zones, so that radiative 
winds and mass loss (boosted by the high surface enrich¬ 
ment in heavy elements) are highly enriched with the CNO 
elements; this cannot be obtained by models of metal-free 
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Table 2. In this Table, we report the numerical values which we predict for the return mass fraction (_R) and the yields of oxygen and 
metals per stellar generation (j/q and 1/Zi respectively) as functions of the metallicity Z and for different IMFs. These values have been 
computed by assuming that the upper mass cutoff of the IMF rriup = 100 Mq. The stellar yields are the ones of Romano et al. (2010). 


z 

R 

yo 

yz 

Stellar yields: Romano et al. (2010) 

R yo yz R yo 

yz 

R 

yo 

yz 


IMF; Salpeter (1955) 

IMF: Chabrier (2003) 

IMF: Kroupa et al. (1993) 

IMF: Kroupa (2001) 

1.0 X 10-1“ 

0.285 

0.028 

0.042 

0.436 

0.059 

0.088 

0.284 

0.015 

0.026 

0.411 

0.053 

0.079 

1.0 X 10-5 

0.285 

0.018 

0.028 

0.436 

0.039 

0.059 

0.284 

0.009 

0.017 

0.411 

0.035 

0.053 

5.0 X 10-5 

0.285 

0.018 

0.028 

0.436 

0.039 

0.059 

0.284 

0.009 

0.017 

0.411 

0.035 

0.053 

1.0 X 10-4 

0.285 

0.018 

0.029 

0.436 

0.039 

0.060 

0.284 

0.009 

0.018 

0.411 

0.035 

0.054 

5.0 X 10-4 

0.286 

0.018 

0.029 

0.437 

0.039 

0.060 

0.285 

0.009 

0.018 

0.412 

0.035 

0.054 

1.0 X 10-5 

0.286 

0.018 

0.029 

0.438 

0.039 

0.060 

0.287 

0.009 

0.017 

0.414 

0.035 

0.054 

5.0 X 10-5 

0.292 

0.018 

0.027 

0.447 

0.038 

0.057 

0.295 

0.009 

0.016 

0.422 

0.034 

0.051 

1.0 X 10-2 

0.295 

0.018 

0.028 

0.451 

0.038 

0.060 

0.299 

0.010 

0.017 

0.425 

0.034 

0.054 

2.0 X 10-2 

0.298 

0.018 

0.031 

0.455 

0.037 

0.065 

0.302 

0.010 

0.018 

0.430 

0.034 

0.059 


non-rotating massive stars (see, for example, Maeder 2009, 
for a detailed discussion). 

On the other hand, yo is strongly dependent on the 
assumed IMF. The highest oxygen yield is obtained when 
adopting a Chabrier (2003) IMF, because this particular 
IMF contains the largest number of massive stars compared 
to the other IMFs explored in this paper (see Fig. 3). The 
IMF of Kroupa et al. (1993), instead, predicts the lowest yo, 
since it contains the lowest fraction of high mass stars. In 
this context, it is important to distinguish the two IMFs 
suggested by Kroupa. In fact the Kroupa (2001) is very sim¬ 
ilar to the IMF of Chabrier (2003) and predicts a substan¬ 
tially higher yield than Kroupa et al. (1993). The Salpeter 
(1955) IMF predicts a net yield roughly halfway between the 
Chabrier (2003) and Kroupa et al. (1993) IMFs. 

In Table 2, we show also how yz (where Z here is the 
sum of all metals) is predicted to vary as a function of the 
different IMFs and metallicities. These values are shown here 
only for reference with previous works attempting to model 
the total metal content of galaxies, however we caution the 
reader against a blind application of analytical models as¬ 
suming the IRA to the total metal content. Finally, in the 
same Table, we report the values of the returned fraction 
R, which is rather constant as a function of metallicity but 
shows some change for different IMFs. The approximate con¬ 
stancy of i? as a function of the metallicity is due to the fact 
that this quantity is strongly dominated by the H and He 
contributions. 

Our results for yo{Z), j/z(Z) and i?(Z), as obtained with 
the Nomoto et al. (2013) set of stellar yields, are reported in 
Table 3. The effect of the various IMFs is the same as dis¬ 
cussed above for the stellar yields of Romano et al. (2010, 
see Table 2). On the other hand, by comparing the predicted 
net yields of metals and oxygen of Romano et al. (2010) with 
the ones of Nomoto et al. (2013), we can quantify the uncer¬ 
tainty introduced by different input stellar yields by a factor 
which is ~ 1.5. 

In Fig. 5, we explore how the choice of the upper cutoff 
of the IMF, rriup, affects the average net yield of oxygen, 
where (yo) stands for the net yield of oxygen as averaged in 
the metallicity range 1.0 x 10“® < Z < 2.0 x 10“^. Our re¬ 
sults are shown for different IMFs (different colors) and dif¬ 
ferent stellar yield compilations (thick and thin lines repre¬ 


sent our results with Romano et al. 2010 and Nomoto et al. 
2013, respectively). By looking at the figure, we see that 
the difference between the case with rriup = 100 Mq and 
rriup = 40 Mq is as much as a factor of about two. The dif¬ 
ference is only ~ 50 per cent between rriup = 100 Mq and 
rriup = 60 Mq. We find that, when assuming the Chabrier 
(2003) and Kroupa et al. (1993) IMFs, the differences in 
{yo) with different upper mass limits are almost doubled 
and halved, respectively, with respect to the case with the 
Salpeter (1955). Fig. 5 shows that the global variation of 
{yo) spanned by all “classical” IMFs and the possible range 
of rriup is nearly a factor of ten. 

By looking at Fig. 5, the curves corresponding to the 
Romano et al. (2010) stellar yields lie always above the 
curves with the Nomoto et al. (2013) stellar yields. This 
difference enlarges as rriup increases, because Romano et al. 
(2010) provide the oxygen stellar yields up to 60 Mq, while 
Nomoto et al. (2013) only up to 40 Mq (see also Fig. 1), and 
we keep these stellar yields constant for stars with larger ini¬ 
tial stellar mass. The latter assumption can introduce a sys¬ 
tematic effect in the final values of yo. We find that, by vary¬ 
ing the upper limit of the integral at the numerator of equa¬ 
tion 2 and by normalizing the IMF up to rriup = 100 Mq, 
the trend of the resulting {yo) is similar to the trend of {yo) 
as a function of the upper cutoff of the IMF (see Fig. 5). 

Assuming a top-heavy IMF can cause an even larger 
increase of the yield of oxygen per stellar generation, being 
larger the number of massive stars which are present. We 
explore the effect of two top-heavy IMFs (Cappellari et al. 
2012; Arimoto & Yoshii 1987), both defined as a single-slope 
power law. Our results for the {yo) vs. rriup relations are 
shown in Fig. 6 for different IMFs and stellar yield assump¬ 
tions. By looking at Fig. 6, as the slope of the IMF is de¬ 
creased from X = 1.35 (Salpeter 1955) down to a; = 0.95 
(Arimoto & Yoshii 1987) and x = 0.5 (Cappellari et al. 
2012), the IMF becomes top-heavier and the value of yo be¬ 
comes larger and larger; furthermore, the standard deviation 
of {yo) becomes larger as the slope x is decreased, indicat¬ 
ing that {yo) is slightly more influenced by the metallicity- 
dependence of the assumed set of stellar yields. 

The predicted values of yo for single-slope top-heavy 
IMFs are very high and they may either imply that a top- 
heavy IMF star formation mode has only lasted for a short 
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Table 3. In this Table, we report the numerical values which we predict for the return mass fraction (R) and the yields of oxygen and 
metals per stellar generation (j/q and 1/Zi respectively) as functions of the metallicity Z and for different IMFs. These values have been 
computed by assuming that the upper mass cutoff of the IMF rriup = 100 Mq. The stellar yields are the ones of Nomoto et al. (2013). 


z 

R 

yo 

yz 

Stellar yields: Nomoto et al. (2013) 

R yo yz R yo 

yz 

R 

yo 

yz 


IMF: Salpeter (1955) 

IMF: Chabrier (2003) 

IMF: Kroupa et al. (1993) 

IMF: Kroupa (2001) 

0.0 

0.261 

0.021 

0.043 

0.403 

0.044 

0.087 

0.244 

0.011 

0.024 

0.380 

0.040 

0.079 

1.0 X 10-3 

0.293 

0.018 

0.026 

0.450 

0.038 

0.055 

0.291 

0.009 

0.014 

0.424 

0.034 

0.050 

5.0 X 10-3 

0.300 

0.016 

0.025 

0.459 

0.034 

0.052 

0.300 

0.008 

0.013 

0.433 

0.030 

0.047 

1.0 X 10-2 

0.302 

0.015 

0.024 

0.463 

0.032 

0.051 

0.303 

0.008 

0.013 

0.436 

0.029 

0.046 

2.0 X 10-2 

0.305 

0.014 

0.023 

0.466 

0.030 

0.049 

0.307 

0.007 

0.012 

0.439 

0.027 

0.044 

5.0 X 10-2 

0.304 

0.017 

0.023 

0.466 

0.036 

0.049 

0.307 

0.009 

0.012 

0.439 

0.032 

0.044 


interval of the galaxy lifetime and not relevant for the inte¬ 
grated metal production, or that a single power-law is not 
a proper representation of the “top-heavy IMF” and that a 
broken power-law is a more proper description. 

There is increasiirg evidence in the literature that the 
IMF may vary among different types of stellar systems, such 
as spheroids and disk galaxies, as well as faint dwarf galax¬ 
ies (e.g. see Cappellari et al. 2012; Conroy & van Dokkum 
2012; Weidner et al. 2013). Observationally, a Kroupa et al. 
(1993) IMF is favoured in describing the chemical evolution 
of the solar vicinity (see Romano et al. 2010), and the disk 
of spirals similar to the Milky Way, while the Kroupa (2001) 
and Chabrier (2003) IMFs are probably better for describ¬ 
ing the evolution of spheroids such as bulges and ellipti¬ 
cals (see, for example, Chabrier et al. 2014). We have shown 
that the range of commonly adopted IMFs (even neglecting 
the extreme top-heavy IMFs) implies a large variation of 
net yield per stellar generation. This fact could add an ex¬ 
tra systematic to studies attempting to model the observed 
abundances, which should be taken into account by prop¬ 
erly using our compilation of yield for the different classes 
of galaxies. More specifically, if the IMF is not universal, we 
can expect a difference in net yield up to a factor larger than 
three for classical, widely-used IMFs, and even much larger 
for top-heavy IMFs. 


4 CONCLUSIONS AND DISCUSSION 

The blooming of extensive spectroscopic surveys of local and 
distant galaxies have fostered the use of metallicities to con¬ 
strain the star formation history, feedback processes (out¬ 
flows) and gas inflows across the cosmic epoch, by compar¬ 
ing the observations with the expectations of analytical and 
numerical models. One key element in such a comparison is 
the yield per stellar generation, which is often assumed as 
a fixed value. In this article, we have used a numerical code 
of chemical evolution, which includes modern stellar nucle- 
osynthetic yields, to explore the effect of the metallicity and 
IMF on the net yield of oxygen per stellar generation and 
on the return mass fraction; our results have shown that 
the yield can change by large factors. Therefore, if the yield 
associated with the appropriate IMF is not used, this can 
produce inconsistent results and large systematic errors. 

We have provided results for two different sets of stellar 
yields, which are Romano et al. (2010) and Nomoto et al. 


(2013). The former is an inhomogeneous set, since it com¬ 
bines results of hydrostatic burning in rotating massive stars 
(He, C, N, and O from Geneva stellar models) with results 
of explosive burning without stellar rotation (heavier ele¬ 
ments from Kobayashi et al. 2006), which is - in principle 
- an incorrect treatment. On the other hand, Nomoto et al. 
(2013) is a homogeneous set of stellar yields, with the only 
limitation of not including the effect of stellar rotation, ex¬ 
pected to have a strong impact at the very low metallicities. 
We have found that the uncertainty introduced by assuming 
different sets of stellar yields can be quantified by ~ 0.2 dex. 

The effect of assuming different IMFs can cause large 
differences in the net oxygen yield. We have found that the 
Kroupa (2001) and Chabrier (2003) predict the highest oxy¬ 
gen yield, roughly a factor of two higher than for a Salpeter 
(1955) IMF. On the other hand, by assuming a Kroupa et al. 
(1993) IMF, we obtain the smallest net yield, roughly a fac¬ 
tor of two lower than for a Salpeter (1955) IMF. 

The yield per stellar generation also depends signifi¬ 
cantly on the upper mass cutoff of the IMF. The differences 
between the case with rriup = 100 Mq and rriup = 40 Mq are 
of the order of a factor of two with the Salpeter (1955). Since 
the IMF of Chabrier (2003) predicts a larger number of mas¬ 
sive stars, that difference is doubled with this IMF, whereas 
it is halved with the Kroupa et al. (1993) IMF, which pre¬ 
dicts the lowest number of massive stars. If one takes into 
account both the variation with IMF shape and upper stellar 
mass cutoff, the variation of the yield per stellar generation 
can span more than a factor of ten. 

We note that populations of highly enriched galaxies - 
whose metallicities were deemed uncomfortably high - can 
be easily explained by means of a large yield per stellar gen¬ 
eration, as the one associated with commonly used IMFs, 
such as Chabrier (2003) and Kroupa (2001). Similarly, our 
results should warn about a proper use of the so-called ef¬ 
fective yield, yeg = Z/ In i which is observationally de¬ 

rived by inverting equation (1). In particular, the finding of 
J/eff < yz is generally modelled in terms of enriched outflows, 
inflow of pristine gas, or both (e.g. Tremonti et al. 2004; Erb 
2008; Mannucci et al. 2009; Troncoso et al. 2014), whereas 
the finding of i/eff > yz is sometimes used as an indica¬ 
tor of inaccurate metallicity measurements or inappropriate 
metallicity calibrations. We conclude that the deviation of 
2/eff from a “fiducial”, “true” yield also may be partly associ¬ 
ated with IMF being different than assumed. For the same 
reason, high values of the the effective yield (and in partic- 
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ular j/eff > yz) may be indicative of an IMF favouring mas¬ 
sive stars (Chabrier 2003; Kroupa 2001) and/or of an high 
mass cutoff of the IMF itself. By assuming single-slope top- 
heavy IMFs, such as the ones proposed by Arimoto & Yoshii 
(1987) or Cappellari et al. (2012), we find very high values 
for the yields of oxygen per stellar generation, which are 
also slightly influenced by the metallicity-dependence of the 
stellar yields. 

The dependence on metalllcity of the yield is reassur¬ 
ingly small. A significant variation is only found at extremely 
low metallicities (Z = 1.0 x 10“^°). Although the latter re¬ 
sult may be an important aspect to take into account for 
models of primordial galaxies, it strongly relies on the as¬ 
sumed set of stellar yields, which are particularly affected 
by uncertainties at extremely low metallicities. 

We remind the reader that IRA provides correct results 
only for chemical elements restored into the ISM on short 
typical time-scales; oxygen represents the best example of 
such a chemical element, since it is also the most abundant 
metal by mass in the Universe. On the other hand, analyti¬ 
cal models working under IRA assumption fail in following 
the evolution of chemical elements produced by long-lifetime 
sources; examples of such chemical elements are carbon, ni¬ 
trogen and iron. Hence we warn the reader against a blind 
application of IRA for the total mass of metals. In order to 
take Into account the lifetimes of the various stellar produc¬ 
ers in detail, numerical models of chemical evolution should 
be used. 

Our compilation of numerical values of the yield per 
stellar generation, for different IMFs, different upper mass 
cutoffs and different metallicities, will hopefully be useful to 
properly investigate the metallicity in galaxies across cosmic 
epochs, by tackling one of the (generally not acknowledged) 
major sources of uncertainty. 
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